Thermophoresis as persistent random walk 
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In a simple model of a continuous random walk a particle moves in one dimension with the velocity 
fluctuating between +v and —v. If v is associated with the thermal velocity of a Brownian particle 
and allowed to be position dependent, the model accounts readily for the particle's drift along the 
temperature gradient and recovers basic results of the conventional thermophoresis theory. 

PACS numbers: 05.40.-a, 66.10.C-, 82.70.Dd 

Thermophoresis [H, 0, S 0] is the systematic drift of a Brownian particle caused by a temperature gradient VT. 
The steady-state thermophoretic velocity acquired by the particle is given by 

Vt = -DtVT, (1) 

where Dt is generally referred to as thermal diffusion coefficient. A similar phenomenon exists on the molecular level 
in gas mixtures and is known as the Soret effect. In most cases Dt is positive, i.e. the particle moves toward a 
colder region. However, thermophilic behavior, Dt < 0, has also been observed experimentally 0]- Among other 
well-known field-driven transport effects, thermophoresis is perhaps the most subtle: it takes place in a stationary 
state when the pressure is uniform, and therefore the average force on a particle which is fixed in space is strictly 
zero. Interesting in itself, thermophoresis has r nany ap plications, in particular in the context of molecular Brownian 
motors driven by thermal fluctuations [1, H, [l^, ■ 

Despite many efforts, a realistic theoretical model of thermophoresis, valid over a wide region of densities of the host 
fluid, seems to be lacking, and even the physical origin of the effect is still under debate [if. Many authors argue that 
the effect cannot be described by a simple modification of the over-damped diffusion equation in the configurational 
space only (Smoluchowski equation), but requires a kinetic description. Well-known example of such approach is that 
by van Kampen p^ . which takes as a starting point the Kramers equation with position-dependent temperature 
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Here f{x,v,t) is the joint probability density of position and velocity, F is an external force, and 7 is the friction 
coefficient which may depend on x and is assumed to be large. Using the expansion in powers of 7~^, one can obtain 
to the lowest order the corresponding overdampcd Smoluchowski equation for the spatial density f{x,t) in the form 

df d ( d dT\ 

where the mobility /i = 1/7717, the diffusion coefficient D = kT fi, and the thermal diffusion coefficient is 

DT^ — ^^^kfi. (4) 
7717 1 

Therefore, the thermophoretic force Ft on the particle, defined through the relation ^Ft — —Dt^T, equals 

Ft = -fcVT. (5) 

The result ^ often underestimates the value of Dt by up to one order of magnitude [l], and cannot, of course, 
account for thermophilic behavior corresponding to a negative Dt- Another restriction of the approach is that 
the validity of the Kramers equation for the case of nonuniform temperature is postulated rather than proved. More 
"microscopic" approaches lead to the expression for Dt with additional terms involving rather complicated correlation 
functions |14|. [isl fioj . However, conceptually the result Q is important as a demonstration that the expression for 
the current J in the Smoluchowski equation ^ cannot be obtained by simple insertion of the position dependence 
into the conventional drift-diffusion expression J = — Dd/dx)f for an isothermal system. 
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The simplicity of the result (|4]) suggests that it might have a transparent physical interpretation. Yet the 
transition from the Kramers equation ^ in phase space to equation ([3]) for the spatial distribution function 
f{x,t) = J dvf{x,v,t) (the elimination of the fast variable v) is not a trivial step [13, [3, which makes the ori- 
gin of the result ^ somewhat obscure. The same perhaps may be said of the Luttinger's method of fictitious external 
fields tia], which also leads to the result (g]) [lO,!!!. The aim of this Letter is to formulate a minimal qualitative 
model of thermophoresis, which leads to the expression ([4]) elementary and directly. Besides pedagogical merits, such 
model might be useful for numerical modeling of stochastic processes at nonuniform temperature. 

The model is a slightly generalized version of the well-known stochastic process of the continuous persistent random 
walk m, [H, [l^l . In its simplest setting, the process describes a particle moving in one dimension with fixed speed 
V suffering occasionally a complete reversal of direction. Let f+{x,t) and f-{x,t) be the probability density for the 
particle moving to the right and to the left, respectively. Reversals of velocity are Poisson distributed, i.e. occurring 
with a constant rate 1/2t, so that the probability for reversal in a time interval dt is dt/2T. For an infinitesimal time 
step one can write 

f+{x,t + dt) = f+{x~vdt,t) - f+{x ~ V dt,t) dt/2T + f_{x + V dt,t) dt/2T, (6) 
and a similar equation for f^{x,t). The corresponding differential equations read 

df+^^dU^ f+-f- ^ ^ ^ , f+-f- .7^ 
dt dx 2r ' dt dx 2t ' ^ ' 

and lead to the telegrapher's equation for the total density f = f+ + .f~, 



dt'^ T dt dx"^ 
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The same equation holds also for the difference A = /+ — /_ , and therefore for the components /+ and separately. 
Suppose the particle at t = is at the origin with equal probability to be in each of the two velocity states, 

Here the second initial condition follows from the first one and Eqs.([7])- Respectively, the initial conditions for the 
total density / = /+ + /- are 

f{x,Q)^6{x), ^%^ = 0- (10) 

The corresponding solution of the telegrapher's equation is well known ^2^. In the long-time limit t ':S> t , vt x it 
coincides exactly with solution of the the Smoluchowski equation 



f[x,t) « (47rDi)-i/2gxp ^- 
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(11) 
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with, the diffusion coefficient D — tv^. However, unlike the overdamped Smoluchowski equation ((S]), equations ([7]) 
incorporates effects of inertia of the particle. This advantage, which was recognized and used beneficially in many 
previous works (see [2^ and references therein), allows to account for thermophoresis in a particularly simple way. 

To link the model to the problem of thermophoresis, it is natural to identify v with a typical thermal speed of 
a Brownian particle. The choice is ambiguous. For instance, one can set v equal to the root mean square velocity 
Vrms = (v^)^^^ = \/3fcr/m. However, one can show that in this case the final result for the thermophoretic force 
would differ from Eq.([S]) by the factor 3/2 (see Eg. p^ below). As will be shown, a perfect agreement with the 
standard results Dt = D/T and Ft = —k'VT can be achieved if v is identified not with Vrsm but with the most 
probable speed of a Brownian particle Vmp (for which the Maxwell speed distribution has a maximum): 



V = Vmp = \/2fcr/m. (12) 

The second parameter of the model l/r should be identified with the friction coefficient 7 which appears in the 
Kramers equation ([2]) and in the corresponding Langevin equation v = —jv + 

For nonuniform temperature, the velocity in Eqs. (jT]) is position dependent, v{x) = ^2kT{x) jm. In this case 
instead of the telegrapher's equation ([8]) one obtains [25| 

92/ 1 a/ f dv\ df 

^'t^+Nx: it- (13) 
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Suppose the temperature gradient is constant, so that T{x) = T + xWT and 

VT 



v{x) = y/2kT{x)/ni = v[l + —X ] , (14) 



where v is the thermal velocity corresponding to the temperature T, v — ^IkT jm. Then the equation ()13|) reads as 



5*2 ' T dt^"" t"") dx-^ +"27 

For a small gradient {yT/T)x <C 1, the equation is simplified to the form 



where the thermophoretic force Ft coincides (thanks to the setting v — \j2kT jm) with the expression ([5]) of the 
standard theory: 

1 VT 

Ft ^ --mv^ — = -kVT. (16) 



It can be shown that in the long-time limit the solution of the equation p5|l with the boundary conditions (fTO|) 
coincides with the solution of the overdamped equation ([3]) and describes the drift along the temperature gradient, 
superimposed on the diffusion 

/(.,,) ^^.xpf-fc^) ,17) 



V 4T>t 

with the drift velocity Vt — Ftt/iti. Using (fTB|) and recalling D — rv^ ^ one gets Vt ~ —Dt^T with Dt = D/T, 
thus recovering the result (|4|) of the conventional theory. 

The asymptotic solution p7p can be obtained as follows. After applying the transformation f{x,t) — 
(j)(x, t) exp(— t/2T + FTx/2mv'^), the equation takes the form 



9*2 ^2 



with l/r2 = l/4r2 — FT/^m^v^, while the boundary conditions corresponding to (fTO|) read 

^(.,0)=^(.), ^ = 

Note that the model makes sense only under the assumption t2 > 0, which guarantees the positiveness of /+ and 
/- [13 ■ The equation ITSl) is the modified telegrapher's equation whose solution is well known ^28*1. Transforming 
back from to /, the result can be written in the following form 

f[x, t) - exp {-^ + p^) (01 +^2 + (19) 



2t 2niv^ 



where the functions (j)i(x^t) are 



01 (a;,*) = ]^[5{x ~ vt) + 5{x + vt)], (20) 
<p2{x,t) = f— ^^(^;^-|x|), 

4WT \ VT^ J 

Mx,t) = ^ . !^ iJJ-^y2t2_xA e{vt-\x\), 

and d(x) is the unit step function. Using the asymptotic form of the modified Bessel functions for large argument 
Ia{x) w -7==e^, it is an easy matter to prove that in the limit of strong damping Ttt/tou <C 1 and of long time 



t ^ T, vt ^ X, the exact solution is reduced to the simple form ([T7]). 
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One interesting problem related to thermophoresis is that of Brownian motors driven by position-dependent tem- 
perature. In particular, the Biittiker-Landauer motor is essentially a Brownian particle diffusing in the periodic 
potential field subject to a spatially inhomogeneous temperature. In this context it is of interest to generalize the 
model to the case of the presence of an external force F. As discussed in ^2^, in this case the equations ([7]) for f±{x, t) 
should be extended as follows 

df± ^ df± f+-f^ ^ (f ^ f \ 

Respectively, the equation (|15p for the total density f = ,f+ + /- is generalized to the form 

d'^f 1 df ^d^f Ft df d [F \ 



dt'^ T dt dx^ m dx dx 
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This equation differs from the overdamped equation ([3]) of the standard theory by the presence of the term ftt- This 
term, related to inertial effects, is unimportant in the long time limit, but may be responsible for wave-like behavior 
at short times j23 |. 

Summarizing, in this Letter we discussed a stochastic process which underlies thermophoresis in a way similar to 
that as the discrete random walk underlies isothermal diffusion. The model leads to the telegrapher's equation whose 
asymptotic solution coincides with the solution of the overdamped Smoluchowski equation ([3]). This is consistent with 
an observation that the telegrapher's equation is reduced to the Smoluchowski equation under conditions u — s- oo, 
T — > 0, v'^T = const. Since the model takes into account inertial effects, one might hope that it can resolve difficulties 
of the overdamped theory of thermally driven Brownian motors [lOl, ITll.ll2.]. However, similar to the Kramers equation 
(12) , the model is based on the assumption that the thermalization of the particle to a local temperature is faster than 
any other process involved. Since the characteristic time scales for inertial and thermalization effects are typically 
the same, the application of the model beyond the overdamped regime, strictly speaking, is not justified, and should 
be undertaken with caution. Despite of this limitation, the mapping of thermophoresis onto a random walk problem 
may offer some advantage for both analytical and numerical modeling, in particular for problems with complicated 
or velocity-dependent ^Q, [IS] boundary conditions. 
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